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ABSTRACT 

Context. XTE J1818-245 is an X-ray nova that experienced an outburst in 2005, as first seen by the RXTE satellite. The source was 
observed simultaneously at various wavelengths up to soft y-rays with the INTEGRAL satellite, from 2005 February to September, 
during our INTEGRAL Target of Opportunity program dedicated to new X-ray novae and during Galactic Bulge observations. 
Aims. X-ray novae are extreme systems that often harbor a black hole, and are known to emit throughout the electromagnetic spectrum 
when in outburst. The goals of our programme are to understand the physical processes close to the black hole and to study the possible 
connection with the jets that are observed in the radio. 

Methods. We analysed radio, (N)IR, optical. X-ray and soft y-ray observations. We constructed simultaneous broad-band X-ray 
spectra covering a major part of the outburst, which we fitted with physical models. Analyzing both the light curves in various energy 
ranges and the hardness-intensity diagram enabled us to study the long-term behaviour of the source. 

Results. Spectral parameters were typical of the Soft Intermediate States and the High Soft States of a black hole candidate. The 
source showed relatively small spectral variations in X-rays with considerable flux variation in radio. Spectral studies showed that 
the accretion disc cooled down from 0.64 to 0.27 keV in ~100 days and that the total flux decreased while the relative flux of the 
hot medium increased. Radio emission was detected several times, and, interestingly, five days after entering the HSS. Modeling the 
spectral energy distribution from the radio to the soft y-rays reveals that the radio flares arise from several ejection events. 
Conclusions. XTE J1818-245 probably belongs to the class of low-mass X-ray binaries and is likely a black hole candidate transient 
source that might be closer than the Galactic Bulge. The results from the data analysis trace the physical changes that took place in the 
system (disc, jet/corona) at a maximum bolometric luminosity of 0.4—0.9 xlO^** erg s"' (assuming a distance between 2.8-4.3 kpc) 
and they are discussed within the context of disc and jet models. 

Key words, black hole physics - stars: individual: XTE J1818-245 - gamma rays: observations - X-rays: binaries - infrared: general 
- radio continuum: general 



1. Introduction 

X-ray Novae (XNe), also known as Soft X-ray Transients, are 
accreting Low-Mass X-ray Binaries (LMXB) that spend most of 
their time in a faint, quiescent state. They undergo large ampli- 
tude outbursts with rise times of only a few days or weeks, with 
typica l recurrence periods of many years dTanaka & Shibazakil 
Il996l) . The picture commonly accepted for an XN involves the 
transition from a low-mass accretion rate state to a high-mass 
accretion rate optically thick accretion flow, triggered by an ac- 
cretion disc instability. The optically thick and geometrically 
thin accretion disc has a varying inner radius and temperature, 
emitting at typical X-ray energies of ~1 keV. This region is 
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probably surrounded by a hot corona, where soft X-ray pho- 
tons originating in the disc undergo inverse Comptonization, 
emitting a power law spectrum up to y-rays. A relativistic jet 
might be present, typically observed in the radio. These spec- 
tral characteristics are coupled to different levels of variability, 
to Quasi Periodic Oscil lations (QPOs) observed in the po wer 
spectrum spectrum (e.g.. lBelloni et al.ll2001l;lBellonill2005al) and 
to changes in the radio. Depending on the relative strengths of 
each component and on how they vary, several spectral state s 
have been identified: see, e.g., iMcClintock & RemillardI ( l2006l) . 
iHoman & Bellonil (|2005|) . In the second classification, the two 
main spectral states are the Low/Hard State (LHS), dominated 
by non-thermal emission and the High Soft State (HSS), domi- 
nated by emission from the accretion disc. In the LHS, the fast 
time variability is dominated by strong (~30% fractional rms) 
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Fig. 1. RXTE/ASM light curve of XTE J1818-245 during the 2005 outburst. The times of the observations are indicated for each 
instrument as well as the INTEGRAL revolutions (lasting ~3 days between consecutive vertical dashed lines). Other observations 
lasting from 100 to 3000 s are indicated (see Table[T]). 



band-limited noise. At times, low-frequency QPOs are present. 
In this state, flat-spectrum radio emission is observed, associated 
with compact jet ejection dCorbel et al.ll2000l 120031: iGallo et alj 
I200I l2006t iFender et alj|2004ft . In the HSS, only weak power 
law noise is present i n the power spectrum. No core radio emis- 
sion is detected (see iFender et al.lll999l 120051 : "quenching" of 
the jet). Further states have been identified as "Intermediate" 
based on the above-mentioned differences in the soft/hard X- 
ray components, variability and radio emission: the HIMS (Hard 
InterMediate State) and SIMS (Soft InterMediate State). In the 
HIMS, the energy spectrum is softer than in the LHS, with evi- 
dence for a soft thermal disc component. The power spectra fea- 
ture band-limited noise with characteristic frequency higher than 
the LHS and usually a rather strong 0.1-15 Hz type-C QPO (see, 
e.g., Casella et al. 2005). In the SIMS, the disc component domi- 
nates the flux. No strong band-limited noise is observed but tran- 
sient type-A and type-B QPOs are seen (the frequency of which 
spans only a limited range). As in the HSS, no core radio emis- 
sion is detected. 

XTE J18 18-245 was discovered by the All-Sky Monitor 
(ASM) telescope on board the RXTE satellite on 2005 August 
12. The ASM Hardness Ratio (hereafter, HR) indicated a very 
soft spectrum, often associated with a B lack Hole (BH) a s 
the compact object of the binary system (iLevine et al.l l2005l) . 
Follow-up i?XTE/Proportional Counter Array (PCA) observa- 
tions provided a refined p osition and no pulsations were detected 
(iMarkwardt et all l2005l) . On August 16-17, the INTEGRAL 
Soft Gamma-Ray Imager (IBIS/ISGRI) observed th e source for 
12.6 k s during the Galactic Bulge (GB) program: IShaw et al] 
(I2005h found a position in the ISGRI mosaic image consistent 
with (and with smaller uncertainties) than the one found by 
the RXTE instruments. Soon a fter, the optical counterpart was 
identi fied dSteeghs et al.l l2005l) and the ^wZ/f/X-ray Telescope, 
XRT, ( Still et al. 20051) improved the X-ray position. Finally, 
iRupen et al.l (l2005h detected XTE J18 18-245 with the Very 



Large Array (VLA) at the optical position: the flux densities at 
4.9 GHz increased from ~7 to 27 mJy between August 20-21. 

Based on the above-mentioned characteristics, 
XTE J18 18-245 was suggested to be a BH Candidate 
(BHC). Target of Opportunity (ToO) observations for XNe in 
outburst were triggered with INTEGRAL, associated with a 
large multi-wavelength campaign. We report here the results 
of the INTEGRAL observations of XTE J1818-245 together 
with Swift, RXTE and NIR/optical/radio data. We start with a 
description of the available data and of the analysis procedures 
employed in Sect. |2l Results are presented in Sect. [3] followed 
by the interpretations and discussions in Sect. ID We summarize 
our conclusions in Sect.|5] 



2. Observations and data reduction 

Table [T| summarizes the observations analysed in this paper, 
giving for each of them the instrument data available, energy 
ranges, dates, exposure times and modes. Figure [1] shows the 
RXTE/ASM soft X-ray light curve of the outburst, with the time 
of the multiwavelength observations indicated. 



2.1. INTEGRAL 

Target of Opportunity observations were performed on (2005) 
August 20-21, 23 and 28-29, corresponding to revolutions 
(hereafter. Rev.) 348, 349 and 35 1 , respectively. We also used the 
data from public observations performed during the GB monitor- 
ing program. 

The IBIS/ISGRI and JEM-X data were reduced with the 
standard analysis procedures of the Off-Line Scientific Analysis 
OSA 7 ■ Q released by IS DC, with algorithms de scribed in 
iGoldwurmetaP (l2003h and lWestergaard et al.l (l2003l) for ISGRI 
and JEM-X respectively. Systematic errors of 2% were added 
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Observatory Instrument Bandpass Period Total Exposure Number Observation Type 









(MJD-53596.5r 


(ks) 


of Obs. 




INTEGRAL 


IBIS/ISGRl 


18-200 keV 


2.5-23.8 


315 


3/4 


ToO (5x5'')/Galactic Bulge" (hex") 




JEM-X 


5-25 keV 


5.9-14.8 


240 


3 


ToO 


RXTE 


PCA 


3-25 keV 


3.1-103.3 


146 


54 


Public 




HEXTE 


20-150 keV 


3.1-103.3 


146 


54 


PubUc 


Swift' 


XRT 


0.3-10 keV 


5.8 


0.13 


1 


Public 


REM' 


ROSS 


550-800 nm 


22.7-76.5 


3.2 


10 


ToO 


NTT"^ 


SUSI-2 


320-900 nm 


9.6 


0.6 


1 


ToO 


VLA= 


L,C,X,U 


1.4-15 GHz 


5.6-196.01 


4.8/15.4 


1/13 


ToO/Public 


VLBA'' 


s/x 


2.3-8.4 GHz 


8.54-8.67 


10.8 


1 


ToO 



Notes: 

a) MJD 53596.5 corresponds to the maximum flux of the source measured during the period covered by RXTE/ ASM. 

b) 5x5 dither pattern around the nominal target l ocation. 

c) Monitoring program of E. Kuulkers et al. (see lKuulkers et alj2007b . 

d) Hexagonal pattern around the nominal target location. 

e) Snap-shot observations with the specified instruments or receivers. 




Fig. 2. INTEGRAL 20-40 keV IBIS/ISGRI reconstructed sky 
image (logarithmic black-body colour scale) of the region 
around XTE J1818-245 during our first ToO (Rev. 348, 70 ks 
exposure). The source appears at a significance level of 18cr 
over the background, the type I X-ray burster Ginga 1826-24 
at 129cr, the microquasar GRS 1758-258 at 86cr and the neu- 
tron star GX 5-1 at 65cr. 



for both JEM-X (in the 5-25 keV range) and ISGRI (in the 18- 
200 keV range). We used the maps, the response matrices and the 
off-axis and background corrections from OSA 7 . 0. For ISGRI, 
we rebinned the standard spectra to obtain between 4 and 10 
spectral points, depending on the brightness of the source. We 
checked that the spectral index did not change by more than 2% 
during a single revolution {INTEGRAL orbit around the Earth 
lasting ~3 days): this allowed us to sum the flux per INTEGRAL 
Rev. to improve the signal-to-noise ratio for both light curves 
and spectra. The ISGRI image and light curves are respectively 
shown in Figs.|2]and[3}a. 

We could only perform JEM-X analysis on our ToO data 
as the source was outside the Field Of View (FOV) of JEM-X 
in the GB monitoring program. Since the RXTE/PCA instru- 
ment observes the source more frequently and has a higher sen- 
sitivity than JEM-X, the JEM-X data were not included in the 
broad-band spectra so as to be consistent over all our data sets. 




b5b80 bJeOO bj520 bJ640 bj66D 53630 5J700 

Fig.3. Light curves of XTE J1818-245 at several wavelengths, 
a: ISGRI light curve (INTEGRAL revolutions indicated above); 
b, c: RXTE/ASM and PCA light curves (flux obtained from 
spectral modeling, see Sect. 13.4b : d: REM and NTT light curves 
in the R filter (magnitudes de-reddened). Dot-dashed lines cor- 
respond to the radio peaks at 4.9 GHz, then at 1.4 GHz (during 
our coverage) while dashed lines locate the spectral transitions 
from SIMS to HSS and then back to the SIMS. 



However, we verified that the best-fit spectral parameters using 
JEM-X and PCA were consistent within the error bars. 



2.2. RXTE 

We analysed all available observations, taken about once every 
two days from MJD~53599 to 53700. Each observation lasted 
between 1 and 3.3 ks. The RXTE data were reduced with the 
HEASOFT software package v6.4, following standard procedures 
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(see lRodriguez et aLll2003h . Energy spectra were only extracted 
from detector 2 of the PCA top layer. 

We used the latest calibration files provided by the RXTE 
GOF (Guest Observer Facility), which include important cor- 
rections applied to the PCA background since 2007 September 
18. To determine the level of systematic errors, we used Crab 
spectra which were fitted with a model consisting of an absorbed 
power law. The photon index was left free in a narrow band. 
Without adding systematic errors in the spectra, the reduced 
chi-squared (hereafter, was well above 1.0. When including 
a 0.8% systematic error, the;^'^^^ dropped to 1.0; this level was 
adopted for XTE J18 18-245 to account fo r uncertainties in th e 
PCA response. This value is consistent with lJahoda et al.l ( l2006h . 

Data from the High-Energy Timing Experiment 
( HEXTE) were reduced in the standard way as described 
in [Rodriguez et al.l ( l2003h apart from the newer version of the 
reduction software. Due to problems brought about by the 
rocking motion of HEXTE Cluster A, the extraction of spectra 
was restricted to Cluster B. Furthermore, due to poor statistics 
in the HEXTE data points in most of the observations, all 
channels were rebinned by a factor of 4. The resultant RXTE 
spectra of a single observation were fitted simultaneously 
between 3-25 keV for PCA and 20-150 keV for HEXTE. Note 
that we also added the IBIS/ISGRI data when available; one 
INTEGRAL revolution corresponds typically to three distinct 
RXTE observations. 

2.3. Swift/XRT 

Afte r the observation that announced a refined X-ray Swift posi- 
tion et al.l2005h . an additional Swift observation took place 
on 2005 August 20. The source was observed for 129 s in pho- 
ton counting mode in order to assess its position. The XRT data 
were processed with the xrtpipeline (vO.9.9) task applying 
standard calibration, filtering and screening criteria. An on board 
event threshold of ~0.2 keV was applied to the central pixel of 
each event, which has been proven to reduce most of the back- 
ground due to either the bright Earth limb or the CCD dark cur- 
rent (which depends on the CCD temperature). For our analysis, 
XRT grades in the 0-12 range were selected. Since the source 
was extremely bright with more than 50 counts s"', the data suf- 
fered from heavy pile-up. In order to overcome this problem, we 
extracted photons from an annular region with inner and outer 
radius of 25 and 80 pixels respectively. 3831 photons were ex- 
tracted in the 0.3-10 keV energy band and used in the spectral 
analysis. 

2.4. Optical and NIR 

Optical and NIR observation s were p erforme d with the Rapid 
Eye Mount (REM) telesc ope jZerbi et al...2001i:IChincarini et al.l 
120031; ICovino et al.ll2004 equipped with the ROSS optical spec- 
trograph/imager and the REMIR NIR camera. Observations 
of XTE J18 18-245 were carried out during 2005 September- 
October {R band), and on October 29 (JHK bands); late time 
observations in the R band were also performed during quies- 
cence on 2007 July 18. Image reduction was carried out by fol- 
lowing the standard procedures: subtraction of an averaged bias 
frame and division by a normalized flat frame. Astrometry was 
performed using the USNOBl.cQ and the 2MAS^| catalogues. 

' http ://w w w. nof s . navy, mil/data/f chpix/ 
^ http://www.ipac.caltech.edu/2mass/ 



Aperture photometry pe rformed with the SExtractor package 
(iBertin & Arnouts|[T996l) for all the objects in the field. The cal- 
ibration was done against Landolt standard stars for the optical 
filters and against the 2MASS catalog for the NIR filters. In or- 
der to minimize any systematic effect, we performed differential 
photometry with respect to a selection of local isolated and non- 
saturated standard stars. Due to the low galactic latitude (-4°) 
and relatively low absorption in the line of sight, doing this anal- 
ysis presented a number of challenges in obtaining the absolute 
flux calibration (difficulties in finding a source-free background 
annulus). The counterpart to XTE J1818-245 was sufficiently 
isolated from neighboring stars and we verified that the corrected 
light curves of the comparison field stars were sufficiently stable 
(within the errors of ~0.05 mag). 

The source was also observed in the optical and NIR with 
the imager SUSI-2 installed on the NTT (New Technology 
Telescope) at La Silla Observatory (private communication: S. 
Chaty). In this paper, we include only the U, B, V, R and / 
photometry to build the Spectral Energy Distributions (hereafter 
SED) simultaneously with the radio/X-ray/y-ray data. 

2.5. Radio 

2.5.1. VLA 

We observed XTE J18 18-245 with the National Radio 
Astronomy Observatory (NRAO) VLA at 1.4, 4.9, 8.4 and 
15 GHz from 23:34 UT of August 22 to 01:35 UT of August 
23 with the VLA in its C configuration. These observations 
were thus simultaneous with our INTEGRAL run. The receiver 
setup included two Intermediate Frequency (IF) channel pairs of 
50 MHz bandwidth each. We used the phase reference calibra- 
tors J181 1-2055 at 1.4 GHz, J1820-2528 at 4.9 and 8.4 GHz, 
and J1751-2524 at 15 GHz. The flux density calibrator was 
J13314-3030 (3C 286). Snapshots of 10 min were obtained at 
1.4, 4.9 and 15 GHz. We also acquired a ~1 hour light curve at 
8.4 GHz. 

Archival data obtained with the VLA in the same configu- 
ration were also retrieved from the NRAO database. A log of 
the observations, the observed fluxes and their associated errors 
can be found in Table |2] All the data were reduced using stan- 
dard procedures within the NRAO aips software package. We 
restricted the 1 .4 GHz data analysis to baselines above 3 kA to 
avoid the extended galactic diffuse emission as much as possible. 
Images with natural weighting were produced and flux density 
measurements were obtained with the aips task jmfit: the abso- 
lute flux calibration is expected to be accurate to the 3% level. 

2.5.2. VLBA 

From 01:00 to 04:00 UT of August 23, we also observed 
XTE J18 18-245 with the NRAO Very Long BaseHne Array 
(VLBA) simultaneously at 2.3 and 8.4 GHz. These obser- 
vations were simultaneous with our INTEGRAL run and 
overlapped with the end of the VLA observations. They 
were conducted using the phase-referencing technique, switch- 
ing between the phase reference calibrator J 1820-2528 and 
XTE J1818-245(separatedby 1°), with cycling times of 2.5 min- 
utes (100 s for the source and 50 s for the calibrator), compatible 
with the expected coherence times. Scans of 2 minutes were ac- 
quired on the fringe finders J1733-1304 and J2000-1748. The 
data were recorded with 2-bit sampling at 256 Mbps at both cir- 
cular polarizations. A total bandwidth of 64 MHz was provided 
by 8 sub-bands. Half of the bandwidth was used at 2.3 GHz 
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Table 2. VLA observation log for XTE Jl 81 8-245 indicating the calendar dates, MJD and INTEGRAL revolutions during which 
the source was observed, together with flux densities at four radio frequencies, the radio spectral index a, errors and the spectral 
X-ray state. 



Date 


Time 


Rev. 


S 1.4 GHz 


S 4.9 GHz 


S 8.4 GHz 


15 GHz 


Spectral index a 


X-ray state 


(yyyy-mm-dd) 


(MJD) 




(mJy) 


(mJy) 


(mJy) 


(mJy) 






2005-08-20 


53602.10 


347 




7.13±0.07 (i, 3.5fr) 








SIMS 


2005-08-21 


53603.19 


348 




27.33±0.06 (d, 7.7cr) 








SIMS 


2005-08-23 


53605.02 


349 


1.5±0.4 


1.49±0.07 


2.3+1.0 (o, 25cr) 


<0.83 


~0 


SIMS 


2005-08-24 


53606.19 


349 






0.60±0.04 






SIMS 


2005-08-27 


53609.19 


350 




1.12±0.07 








SIMS 


2005-08-29 


53611.19 


351 


9.3±0.5 (d, 3.10-) 


3.72±0.08 


2.27±0.07 




-0.8 ±0.1 


SIMS 


2005-09-04 


53617.09 


353 


<1.22 


0.77±0.08 


0.40±0.09 




-1.2 ±0.5 


HSS 


2005-09-11 


53624.11 




<0.86 


<0.23 


<0.17 






HSS 


2005-09-18 


53631.15 






<0.17 








HSS 


2005-09-23 


53636.12 






<0.13 








HSS 


2005-10-06 


53649.02 






<0.22 








HSS 


2005-10-13 


53656.02 






<0.16 








SIMS 


2005-10-17 


53660.84 






<0.23 








SIMS 


2006-02-26 


53792.51 






<0.23 








Unknown 



Note: 

The large errors at 1.4 GHz are due to bright sources in the primary beam (field of view). Upper limits are at the 3(T level. The large errors at 
8.4 GHz during these observations (Rev. 349) reflect the variability of the source (see text). For the runs where variability is present, we quote in 
parentheses whether the source flux density is increasing (i), decreasing (d) or oscillating (o), together with the significance of this variability. 



and the other half at 8.4 GHz. The data were processed at the 
VLBA correlator in Socorro, using an integration time of 2 s. 
Post-processing and data reduction were conducted using stan- 
dard procedures within the NRAO aips software package. 

3. Results 

The following two sub-sections refer to the hard X-ray image 
shown in Fig. |2] and the optical, (N)IR and X/y-ray light curves 
shown in Fig. [3] 

3. 1 . X-ray position 

In the combined IBIS/ISGRI images obtained near the peak of 
the outburst (August 20-2 1), XTE J 1 8 1 8-245 was detected at 1 8 
and 1 2 cr in the 20-40 and 40-80 ke V energy bands, respectively. 
Figure|2]shows the 20-40 keV IBIS/ISGRI mosaic obtained dur- 
ing ~70 ks within Rev. 348. The best-fit position in this image 
is aj2000 = 18^ 18™ 24!36 and Snooo = -24° 32' 34'.'3, with a n 
accuracy of 1™58 at the 90% confidence level dGros et alJl2003h . 

This is compatible within errors with the most precise posi- 
tion obtained from the XRT data using xrtcentroid: ajaooo - 
18^ 18™ 24^;02 and 6mm = -24° 32' 19'.'3, with an accuracy of 
3'.'5 at the 90% confidence level. 

In turn, this lies just 5'.'8 away from the cen t er of the ra- 
dio position at 4.9 GHz reported bv iRupen et al.l (l2005h at the 
beginning of the outburst; a'j2ooo = 1 8'' 18*" 24';43 + 0:2 and 
5j2000 = -24° 32' 17'.'96 + 0'.'4 (see Sect. [33T2l for an updated ra- 
dio position). The high-energy source and the optical/radio coun- 
terparts are therefore all unambiguously associated with the new 
X-ray transient source. 

3.2. X-ray and soft y-ray iight curves 

Figure [3}a shows the 7iVTEGRAL/IBIS light curve of 
XTE J1818-245 during the outburst. The peak seen by ISGRI 
occurs at approximately MJD 53599: the source flux reaches 



the maximum value of ~1 counts s in the 20-200 keV range, 
which corresponds to ~33 mCrab. Note that due to incomplete 
ISGRI coverage at the beginning of the outburst, the actual hard 
X-ray maximum may have occurred a few days earlier Up to 
MJD 53620, the IBIS/ISGRI count rate slowly decreased. 

Figure Ob shows the 1.5-12 keV RXTE/ ASM daily av- 
erage light curve from (2005) July 29 to December 16, with 
a peak flux around MJD 53596.5. Assuming an exponential 
shape for both the rise and the decay phases seen in the ASM 
light curve, we obtained time constants of, respectively, 5+1 
and 19±1 days. An interesting feature is present between MJD 
53620 and 53635 (Fig. [3]-c): the exponential decrease of the 2- 
20 keV flux ceases, and then begins to increase again during the 
subsequent ~10 days. A shorter event, with smaller amplitude, 
might have also occurred between MJD 53603 and 53606. This 
feature, sometimes known as a secondary maximum or a bump 
has already been noticed in several XNe (e.g., A0620-00, Nova 
Muscae, 4U 1543-47) and wifl be discussed Sect.|4TT] 

3.3. Hardness intensity diagrams and quasi-periodic 
oscillations 

To get a first idea of the spectral behaviour of the source, we 
produced a Hardness-Intensity Diagram (HID) with RXTE/PCA 
Fig, in similar to those widely used in the literature 
Fender et al.ll2004l) . While there was no coverage from the start 
of the outburst to the peak, XTE J18 18-245 did not follow the 
usual path of XNe in outburst on this diagram: although it cer- 
tainly went from high flux and relatively low HR to lower flux 
and higher HR, trac ing the end of the usual Q-shape (see, e.g., 
lBelIonietal.ll2005bl) . there were slight deviations. Indeed, the 
source returned along its path twice and traced a Z pattern in 
each spectral state, evolving from low to high HR (see Sect. 14. It . 

Using this HID obtained with the PCA, together with the 
analysis that will follow (Sect. [34l Fig. |5]and|6l), we identified 
two distinct spectral states. The upper branch on the HID had a 
HR ~0.3, a power law flux representing ~ 10-20% of the total 
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Hardness ratio (6.0-19.1 keV / 2.3-6.0 keV) 

Fig. 4. Hardness-intensity diagram of XTE J1818-245 on a logarithmic scale as obtained with the top layer of RXTE/PCA (2.3- 
19.1 keV, detector 2). The source globally evolved from top to bottom and from right to left then right again, but returned along 
the track twice. The filled squares represent individual RXTE observations, the thicker paths indicate the times of the INTEGRAL 
Rev. (approximately indicated). The SIMS and HSS states are delimited by solid lines perpendicular to the source path. The times 
of radio observations are represented by the open triangles. The big circles correspond to the times of the radio flares (top: 4.9 GHz; 
bottom: 1.4 GHz). The downward arrow on top indicates the HID position during our simultaneous radio and ToO observations of 
MJD 53605 (Rev. 349). The upwai'd arrow marks the HID position during the HSS radio detection on MJD 53617. 



2-20 keV flux, and a spectral index in the range 2.2-2.5 (Fig.|7|l. 
This Soft Intermediate State (SIMS) evolved to a pure HSS 
around MJD 53612: the power law fraction dropped to ~5%, 
with a HR around 0.15. After MJD 53654, the source returned 
to the SIMS, and slowly evolved to harder states: the power law 
fraction increased again to 20%, while the HR increased to 1, 
and the spectral index gradually decreased. Unfortunately, the 
X-ray coverage did not extend beyond MJD 53700, and thus 
there is no proof that the source eventually evolved to the LHS 
(the photon index was ~2.3 instead of 1.3-1.5 and the power 
law fraction too l ow compared to the high values -greater than 
80%- reported in iHoman & Bellonil 120051: see Table [3]l, as is 
usually the case for such objec ts (e.g., ICadoUe Bel et al.ll2004t 
iMcClintock & Remillardll2006l) . In this work, the last and hard- 
est data we have analysed occur when the source was still in the 
SIMS. 

The RXTE PDS showed very weak variability, with the 
power decreasing against frequency in a roughly power law-like 
fashion below a few Hz. The total RMS level was around 6-8% 
before MJD 53608 and too faint (below the detector noise) to 
be quantified. Above a few Hz, the PDS was quickly falling to 
the level of the white noise. QPOs were also searched for, but 
none were detected. To quanti fy this, we used the following for- 
mula from Ivan der KlisI (l2006h which reflects the detection limit 
at which a QPO would be picked up at the 3 cr level: 

, L s+B IfwhmY 



where A is the detection level (at 3cr), S the source count 
rate, B the background level, FWHM the width of the searched 
QPO and T the exposure time in seconds. In the case of 
XTE J1818-245, for a typical QPO with a 1 Hz FWHM, the 



limit at which a QPO would have been detected at the 3cr level 
is 2% rms amplitude at the onset of the outburst, and 8% when 
the X-ray flux was minimum (at the end of the outburst). The 
absence of strong QPOs and of strong band-limited noise is ex- 
pected during soft states like the HSS and the SIMS, giving fur- 
ther support to the proposed spectral classification below. 

3.4. X-ray and j-ray spectra 

Using XSPEC vll.3.2 (lArnaudlll996h . we fitted spectra from 
the PCA (3-25 keV), HEXTE (20-IOOkeV), IBIS (18-200 keV) 
and XRT (0.3-10 keV) instruments simultaneously. Several 
models were tested when analyzing the spectra. A normaliza- 
tion constant was added to account for uncertainties in the cross- 
calibration of the instruments. 

The data were well fitted using an absorbed power law 
combined with a multicolour black-body and a fluorescent 
(Gaussian) iron line. Neither a reflection component nor a high- 
energy cut-off" were needed at any time in our fits, at least up 
to 150 keV: this was tested by adding these additional contri- 
butions, and, as we did not observe any improvement in the 
goodness of the fit, we concluded they were not necessary. 
The absorption was fixed to the value found by XRT (A^h = 
5.4^g g X 10^' cm"^) since it cannot be constrained by either the 
PCA or JEM-X data. This value is compatible, within the errors, 
with the average galactic column density in the source direction 
(4.8 X 10^' cm"2, estimated from lSchlegel et al.|[T998h . An iron 
line was necessary to adequately fit the data, but the line cen- 
troids could not be constrained by the PCA and we thus forced 
them to have energies above 6 keV. 

For the disc component, the ezdiskbb model in XSPEC nota- 
tion (IZimmerman et al.l2005l) was used: it assumes a zero-torque 
boundary condition at the inner edge of the accretion disc, which 
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Table 3. Best-fit spectral parameters over all our INTEGRAL observations. 
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Notes: 

Model applied in XSPEC notations: constant*wabs*(ezdiskbb-i-gaussian-i-powerlaw) with A^h fixed to 5.4 x 10^' cm"^. Errors are given at the 90% 
confidence level {At(~ = 2.7). 

a) Disc normalization K is proportional to (R/D)^ cos 9, where R is the inner disc radius in km, D is the distance to the source in kpc and d the 
inclination angle of the disc. 

b) Computed in the 2-20 keV range. 

c) Extrapolated in the 0.001 keV-10 MeV range. 

d) Based on the RXTE/PCA light curve. 
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Fig. 5. Energy spectra of XTE J18 18-245 during the 
first INTEGRAL ToO (Rev. 348, MJD~53602.5) with the 
Swift/XRT (blue), RXTE/PCA (black), RXTE/HEXTE (red) 
and INTEGRAL /IBIS /ISGm (green) data, along with the best- 
fit model: an absorbed multicolour black-body disc and a power 
law with a Gaussian component (see Table [3] for parameter val- 
ues). Residuals (AS x^) in cr units are also shown below. The 
source was in the SIMS. 
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Fig. 6. Energy spectra of XTE J18 18-245 during the third 
INTEGRAL ToO (Rev. 351, MJD~53611.1) with the PCA 
(black), HEXTE (red) and IBIS/ISGRI (green) data, along with 
the best-fit model: an absorbed multicolour black-body disc and 
a power law with a Gaussian component (see Table[3]for param- 
eter values). Residuals (AS x^) in o" units are also shown below. 
The source was in the HSS. 



is not the case in the widely used diskbb model. This condi- 
tion normally leads to a more accurate determination of the in- 
ner radii of the accretion disc, as well as a more physical value 



for the maximum temperature in the disc. The quality of the fit 
to the data was similar for the diskbb and ezdiskbb models: we 
only found a factor ~2.2 of difference between the inner radii. 
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Time (MJD) 

Fig. 7. Spectral characteristics of XTE Jl 8 18-245 during out- 
burst; a: the ASM light curve is reproduced here in order to 
easily compare with the changes in the spectral parameters; b, 
c, d: disc temperature (in keV), photon index and inner radius 
of the accretion disc (in km) derived from the spectral analysis 
(the inner radius values after MJD 53680 were ignored because 
of the high level of uncertainties); e: fraction of the power law 
flux to the total flux in the 2-20 keV range. Dot-dashed lines 
correspond to the radio flares at 4.9 and 1.4 GHz while dashed 
lines indicate the spectral transitions from and to the SIMS (see 
Table[3]and text for details). 



which was consistent with the value found by IZimmerman et aP 
(I2005h . 

The power law component is taken to be a purely phe- 
nomenological model and could signal the presence of a com- 
pact jet, a corona, or reprocessed hard X-ray emissi on due 
to X -ray heating from an extended central source dHvnesI 
2005h. We tried the more p hysi cal Comptonizati on models of 



Sunyaev & TitarchukI (Il980l) and lTitarchud (11994 without sue 
cess (the fit was worse than with the simpler model). We 
also applied much more complex models better adapted for 
such soft states such as compps and eq pair, developed re spec- 
tively by iPoutanen & Svenssoiil ( 1 19961) and ICoppil d 19991) (see 
ICadoUe Bel et al.l (l2006l) for discussions). Unfortunately, the 
source was not bright enough to constrain the spectral parame- 
ters: the resulting did not decrease. Besides, a simple power 
law fit allowed us to compare more easily the spectral parameters 
between the bright beginning and the faint end of the outburst. 

Thus, we fitted all the data with the unique con- 
stant*wabs*(ezdiskbb+gaussian+powerlaw) model (in XSPEC 
terminology) throughout the whole outburst. For all spectra, nor- 
malization constants were in the 0.7-1.4 range for instruments 



other than the PCA (which was frozen to 1). Two examples of 
fitted spectra obtained during our INTEGRAL ToO observations 
are shown in Fig.|5]and|6](in photons cm - s ') with residuals. 

Table |3] lists the spectral parameters derived for all the 
INTEGRAL observations. A plot of the main spectral parameters 
for all the observations is shown in Fig. |7] together with the 
ASM light curve. The spectral parameters were compatible 
with the source being first in the SIMS, then transiting into the 
HSS and returning to the SIMS, this time with a lower flux and 
a higher HR; the dashed lines of Fig. [T] indicate both spectral 
state limits. Throughout the outburst, a regular decrease of the 
disc temperature from a high value of ~ 0.64 to 0.27 keV was 
observed. 

Although not extremely well constrained by the PCA re- 
stricted to energies above 3 keV, the temperature trend could be 
derived: we verified that the disc parameters we derived for the 
XRT data plotted in Fig. |5] are compatible with those epochs 
for which no XRT data is available. The observed decrease for 
the inner disc temperature corresponds to what is usually ob- 
served for a transient source in outburst: the disc dominates the 
emission at the beginning and then gradually cools. The rather 
low temperatures are more consist ent with a binary system con- 
taining a BH than a neutron star (iTanaka & Lewinlll995l) . The 
disc normalization slowly increased during the outburst, from 
~ 150 to 400 xD V( 1 / cos 6) km (D is the distance to the source 
in units of 10 kpc, 9 the inclination angle of the disc). Although 
some caveats may be raised about the exact value of the inner ra- 
dius of the accretion disc obtained from the fits, a clear trend is 
seen: the disc receded globally during the decay phase of the out- 
burst. Again, this is a typical feature in X-ray Novae : during the 
hardening, the disc usually m oves outwards slowly dChen et al.l 
ll997l;ICadolle Bel et al.ll2004l) . even though this is not necessar- 
ily true for all transient sources (see, e.g.. Miller etal. 200a). 
However, this point is still strongly debat ed dDone et al.l 120071 ; 
lRvkoff^etal.1 120071; iGierliriski etakl l2008h as, for example, in 
XTE JI817-330. We do not address this question specifically 
here as our data - with a lower energy limit of 3 keV - are not 
sensitive enough for this purpose. 

At the same time, at higher energies, the power law photon 
index varied between 2.0 and 3.5, thereby confirming that the 
source remained in the softer states. Towards the middle of the 
observations (MJD 53635-53655) the fraction of the power law 
flux compared to the total flux decreased from 5 to 2% (HSS), 
while the photon indices displayed large error bars (statistical 
fluctuations). Subsequently, the power law flux increased rapidly 
from 2 to 5% and gradually returned to the original ~20% level 
(SIMS). The spectral fits and the HR indicate that the spectrum 
was harder at the end of the outburst compared to the peak. This 
is not unusual as XNe tend to evolve towards harder states in 
the late stages of their outbursts. For example, during the decay 
of their outbursts, XTE J1720-318, XTE J1650-500, SWIFT 
J1753. 5-0127 and GX 339-4 showed a slowly receding disc 
with decreasing inner temperature, while at the same time the 
relative amount of the power law contribution increased. 

3.5. Multiwavelength studies 
3.5.1. NIR/optical 

The standard photometry performed with EMMI/NTT on 
XTE J1818-245 around August 24 (between 3-5 UT, i.e.. Rev. 
349) provided the following optical magnitudes (private com- 
munication: S. Chaty): U = 18.24 ± 0.03, B = 18.47 + 0.014, 
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V = 17.46 + 0.01, R = 16.95 + 0.01 and / = 16.3 4 ± 0.01. They 
were p erfectly consistent with the data reported bvlSte eghs et al.l 
(120051) and the REM results which cover a much longer period. 

Figure [3}d shows the magnitude in the R filter cor- 
rected for reddening, obtained with the REM and NTT tele- 
scopes along the outburst. To estimate the Eb-v parame- 
ter necessary for the de-reddening, we used the rel ation 
Nh/Eb-v = 5.8 X 10^1 cm-2 mag ' jBohlin et al.lll978l) and 
the absorption column density estimated from Swift/XKI data 
(A^H=5.4;^gg X 10"' cm"^). The resulting colour excess was 
Eb-v - 0.9 3'^Q''|t mag. Using a standard extinction curve from 
iFitzpatrickId 1999b . we thus obtained the de-reddening parame- 
ters for each of the optic al/NIR filters. 

The source faded consistently, as expected given that the 
REM observations started during the decay phase of the out- 
burst. Some two months after the peak of the outburst (around 
MJD 53672), the magnitudes in the R, J, H and K filters were 
respectively: > 18.3 (at the 3o- Hmit), 16.2+0.3, 14.7+0.2 and 
15.1 ±0.2. This is compatible with the classical picture of XN 
transients in outburst: well after the peak of the outburst, when 
the soft thermal component has faded, the accretion disc still 
contributes (albeit marginally) to the emission in the optical/NIR 
while the jet activity may have ceased. In quiescence, the opti- 
cal/NIR is mainly the quiescent, likely optically thin disc, plus 
the mass donor star, and eventually synchrotro n emission from a 
persis tent jet: see, e.g., SWIFT J1753.5-0127 (ICadolle Bel et al.l 
120071 although this source re mained in the LHS), and more gen- 
erally BH XNe in outbu rst (IChen et al.l 1 19971: iBrocksopp et al.1 



ly^ 

UOOUlHvnes et alj[2002l) as discussed hereafter. 



3.5.2. Radio 

We show in Table |2]the flux densities measured during the VLA 
observations from MJD 53602.10-53792.51. The large errors at 
1 .4 GHz are due to the presence of bright radio sources in the 
field. The radio light curves at 1.4, 4.9 and 8.4 GHz are shown 
together with the ASM X-ray light curve in Fig. [8] 

At least two distinct radio flares were visible (indicated by 
dot-dashed lines in Figs.|3]|2]and[8]i during the period covered by 
our data. First, between MJD 53602 and 53603, XTE J18 18-245 
showed a strong variation at 4.9 GHz: the fl ux density increase d 
from ~1 to 27 mJy in less than 24 hours dRupen et alj 120051) . 
Afterwards, between MJD 53605 and 53611, the flux density 
at 1.4 GHz increased from 1.5 to 9.3 mJy, while at 8.4 GHz it 
only increased from 1.5 to 3.7 mJy. Meanwhile, small spectral 
changes were occurring in the X-ray range : XTE J 1 8 1 8 -245 was 
slowly going from a SIMS to a pure HSS without following the 
classical path in the HID. The X-ray photon index F was ~2.3 at 
that time. 

Regarding the observations conducted on 2005 August 23 
(MJD 53605.02) simultaneously with the INTEGRAL ones 
(Rev. 349), the 8.4 GHz flux density measured in the image 
obtained with the whole data set was 2.37+0.04 mJy. We split 
the data set in 10 and 5 minute intervals, constructed images 
and measured the corresponding flux densities to search for 
variability. The results are shown in Fig. |9] As can be seen, 
the flux density at 8.4 GHz increased from 1 to 4 mJy in just 
30 min and decreased again to the initial value on a similar 
timescale. Meanwhile, the phase-reference source J1820-2528, 
observed every 10 minutes, was clearly stable (with a less than 
7% probability that the source was variable): this showed that 
XTE J1818-245 was variable. The average of the 5-min mea- 
surements is 2.3±1.0 mJy, and this is the value we report in 




53590 53600 
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Time (MJD) 
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Fig. 8. From top to bottom: zoom on the ASM X-ray light curve 
(a) and radio flux densities at three frequencies (1.4, 4.9 and 8.4 
GHz, respectively in panels b, c and d) obtained with the VLA 
from MJD 53590 to 53640. Dot-dashed lines correspond to the 
radio flares at 4.9 and 1.4 GHz, while the dashed lines indicate 
the spectral transition from the SIMS to the HSS and then back 
to SIMS. 



Table |2] Since the source was highly variable and the multi- 
frequency radio data were not strictly simultaneous, it was not 
possible to derive a reliable spectral index for the source. We 
analysed separately the flux densities measured in different IFs 
for the 10 min scans around the maximum at 8.4 GHz, but the 
large errors did not allow us to extract any conclusion on the 
spectral shape. However, comparing the quasi-simultaneous 8.4 
and 4.9 GHz data in Fig. |9] one can see that the spectrum was 
most likely close to flat. 

Similar variability analysis splitting the data of each run 
were conducted for the runs with flux densities above 1 mJy. 
Variability can only be claimed in the following four cases (see 
also the caption of Table |2]i. The source flux density at 4.9 GHz 
was increasing during the first radio observation, with a signifi- 
cance of 3.5cr. It was already decreasing (7.7cr) during the sec- 
ond one, when the first radio outburst took place. It was oscillat- 
ing (25cr) at 8.4 GHz during the third run, as noted in the above 
paragraph. It was marginally decreasing (3. 1 cr) at 1 .4 GHz at the 
time of the second radio outburst. 

The radio position reported by iRupen et al.l ( l2005h was ob- 
tained at 4.9 GHz at the beginning of the outburst: however, at 
that time, the source was highly variable and the flux density de- 
creased from 28.1 +0.1 to 26.7 + 0.2 mJy in just 20 minutes of 
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Fig. 9. Radio flux density measurements as a function of time 
obtained on 2005 August 23 at 8.4 (circles), 4.9 (squares), 1.4 
(diamond) and 15 GHz (lower triangle). The bin size is 5 minutes 
for the first two frequencies and 10 minutes for the other ones. 
The 15 GHz data point corresponds to the 3cr upper limit. 



SWIFT J1753.5-0127 (ICadoUe Bel et all l2007h . This clearly 
supports the hypothesis that XTE Jl 8 18-245 is a good BHC. 

Interestingly, secondary maxima or bumps were seen in the 
ASM and PCA light curves: a smaller one around MJD 53605 
(just before the first observed radio outburst when the disc tem- 
perature was at a relative minimum) and a clear one around 
MJD 53625 (Fig.|7]i. During those bumps, the X-ray flux levelled 
off and was accompanied by a simultaneous increase in the disc 
temperature (Fig. |7]i. Several mechanisms have been proposed 
to explain similar features, one of the most pr omising being 
probab ly the outer disc irradiation proposed by iKing & Rittej 
(Il998h . In this model, the central X-ray source irradiates the ac- 
cretion disc, and maintains it hot until the central accretion rate 
is sharply reduced. Thus, part of the disc keeps a high effective 
viscosity during the outburst. If the irradiation is strong enough 
to heat the whole disc, the X-ray light curve will decay expo- 
nentially. Moreover, in this picture, the secondary maxima are 
explained by a more marked increase in viscosity in the outer 
regions than further in, causing a pulse of extra mass to move 
inwards. Depending on the radial size of the disc region that is 
involved, the secondary maximum thus appears roughly on a vis- 
cous time-scale after the initial o utburst. 

Using the model developed in lShahbaz et al.l (Il998l) . one can 
try to estimate the distance to the source with the following for- 
mula: 



observations. Therefore, to obtain the best position of the radio 
source, we selected the VLA 8.4 GHz data set acquired when 
the source was brighter and stable, which corresponds to 2005 
August 29 (MJD 53611.19). The best-fit position was q'j200() - 
18M8™ 24H30±0^;003 and ^jzooo = -24° 32' 17'.'91 ±a.'07. This 
position is more precis e than (and compatib le within the errors 
to) the one reported by iRupen et alJ (l2005h . The offset is 0'.'42 
only fr om the R-band optical position reported bv lSteeghs et alj 
(l2005b and 0'.'49 from our REM /J-band results, both given 
with an estimated uncertainty around, respectively, (y.'3 and 0'.'2. 
Therefore, these optical positions are fully compatible, within 
the errors, with our radio results. The radio source was unre- 
solved at all frequencies and epochs. 

The VLBA observations conducted at the end of our simul- 
taneous VLA run show that no source was detected within 1" 
of the VLA position of XTE J18 18-245. We could set an upper 
limit of 1.0 mJy beam ' at both 2.3 and 8.4 GHz frequencies. 
The beam at 2.3 GHz had 1 1.6 x 5.9" in position angle (north to 
east) of -5°while the beam at 8.4 GHz had 2.5 x 1.5" in position 
angle of -3°. 



4. Discussion 



Dkpc = 4.3 X 10- 



4. 1 . X-rays and j-rays 



In Sect. 13.31 and 13.41 we showed that, during the broad-band 
observations of XTE J18 18-245 presented in this work, the 
source switched between the SIMS and the HSS, as shown 
by the best-fit parameters obtained in Table [3] and summa- 
rized in Fig. |2] In the X-ray range, the high peak luminos- 
ity, the fast rise timescale, the bright state and the spectral 
properties are typically observed in other dynamically con- 
firmed BHs during their outbursts. The characteristic decay 
time derived is compatible with the us u al behaviour of X Ne 
in outburst jTanaka & Shibazakil [19961 IChen et al.l Il997h as 
seen recently in XTE J 1720-3 18 dcadolle Bel et alJl2004 and 



7] is the radiation efficiency parameter (typically 0.10 for BH 
systems); / is the ratio of the disc mass at the start of the out- 
burst to the maximum possible (we took 0.5 and 1); is the time 
in days when the secondary maximum occurs after the outburst 
peak; Fp is the peak flux in the 0.4-10 keV band, and Td is the 
e-folding time of the decay in days. For XTE J1818-245, we 
obtained - 18-20 d, f , a; 30 d and a distance between 2.8^.3 
kpc. We used it in the luminosity estimation, since this is one of 
the few constraints we have on the source distance which relies 
on a model. This model also leads to the determination of the 
outer radius of the disc (2.5-^.1 xlO'° cm) and its outer viscos- 
ity (0.8-2.3 xlO''' cm^ s"'). It should be noted, however, that 
this model has been tested on ~ 10 sources only, and should not 
be considered as a precise measurement of these parameters. 

The evolution of the source in the HID (Fig. |4]i did not 
strictly follow the usual Q-path of BH XNe: between its first time 
in the SIMS and HSS, the source returned twice along the usual 
path towards higher HR, tracing a small Z-shape, before going 
back to a smoother evolution in the last part of the SIMS. This 
behaviour was first seen in XTE J 1550-564 by iHoman et al.l 
(120011) . and that led these authors to suggest the need for an ad- 
ditional parameter to the accretion rate to drive state changes. 

While our spectral information usually started above 3 keV, 
leading to a possible underestimation of the bolometric lu- 
minosity, we derived for the spectra obtained during the first 
INTEGRAL ToO with Swift/XRT a high unabsorbed 0.3- 
10 keV flux of 3.8x10-** erg cm^^ s-\ The 10-200 keV flux 
reached 2.6x10"^ erg cm"^ s"' while the bolometric flux (ex- 
trapolated from 0.001 keV to 10 MeV) was 4.2x10"^ erg cm^^ 
s"' (dominated by the disc emission at the beginning of the out- 
burst). Assuming a distance to the source between 2.8^.3 kpc, 
this corresponds to an unabsorbed bolometric luminosity of 0.4- 
0.9x10"'^ erg s"' , below the Eddington regime for a 3 solar mass 
BH (i.e., 3.9x10"'^ erg s"'). This high luminosity is not surpris- 
ing considering the soft spectral state at this time. The disc and 
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Fig. 10. Unabsorbed spectral energy distribution of 
XTE J18 18-245 from radio (VLA at 4.9 GHz) to X-ray/soft 
7-ray data (XRT, PCA, HEXTE and ISGRI from 0.3 up to 
200 keV) during the end of Rev. 347/beginning of 348 (i.e., first 
INTEGRAL ToO). A strong increase in the VLA flux density 
was observed in less than 1 day. 



bolometric fluxes then decreased while the high-energy emis- 
sion, well modeled by a power law, increased up to the fraction 
of 25% (see Tableland Fig.U}- 

In the latest stages of the observations presented in this 
work, the source flux was still higher than the expected flux in 
the LHS: the source stayed far from the "hard section" zone 
of the HID. The relatively high power law index of ~ 2.3, 
its low flux contribution and the HR value in the HID clearly 
show the source had not reached the LHS yet; extrapolating 
the later behaviour in Fig. |4] shows that it could have taken 
weeks before the source did so. Until the end of our observa- 
tions, the luminosity stayed greater than the usual quiescence 



values (between 10^°- 10-^^ erg 



XTE J1550-564 



e.g.: 

and references therein); 



er g s i n, 

and H 1743-322: ICorbel et alj|2006l ai 
later, the source did eventually return to quiescence. Several 
studies on BH LMXBs have shown that such sources can 
spend weeks to month s or even years in a giv en spectral state 
(e.g., XTE J1720- 3 1 8. ICadolle BefetaLHIoSl; XTE J65 0-500. 
ICorbel et alj 12004', SWIFT J1753. 5-0127. ISoleri et all |2008[) . 
The spectral evolution of XTE J 1 8 1 8-245 is therefore quite typ- 
ical when compared to other BH LMXBs. 

4.2. Spectral Energy Distribution (SED) 

Two SEDs composed of simultaneous radio, optical, NIR and 
X-ray/soft y-ray data are shown in Figs. [TOl and [TTI In Fig. [TOl 
we plot the data obtained during the first INTEGRAL ToO (Rev. 
348): strong variation in the VLA flux at 4.9 GHz was observed, 
with the flux increasing by ~20 mJy in less than one day. For 
the second SED (Rev. 349, Fig. fTTl i. three radio data points are 
shown at 1.4, 4.9 and 8.4 GHz (the 15 GHz upper limit was 
omitted for clarity) and we added the NTT magnitudes corrected 
for absorption (as described in Sect. 3.5.1). 

Firstly, when fitting separately the radio and the NIR/optical 
data with a simple power law, we obtain clearly distinct slopes. 
From the SED obtained with the most complete coverage (Rev. 
349 , Fig. [TTI ). we derived a radio spectral index a (where 
Sy cc V') close to for the (non-strictly simultaneous) multi- 



Fig. 11. Unabsorbed spectral energy distribution of 
XTE Jl 8 18-245 including radio (VLA at 1.4, 4.9 and 8.4 GHz), 
NIR/optical (NTT) and X-ray/soft y-ray data (PCA, HEXTE 
and ISGRI instruments from 3 up to 200 keV) during the second 
INTEGRAL ToO (Rev. 349). A simple power law is plotted 
through the radio data to show the optical excess: at least 3 
distinct contributions (jet, disc, hot medium) are needed to fit 
the data. 



frequency radio data, while the NIR/optical spectral index was 
-hI.O. In the SED plot the spectral photon indexes F (where 
vFy oc v""^) were thus close to -1 and -2, respectively. Extending 
the power law spectrum derived from the radio data points up to 
the NIR/optical frequencies with the same slope, we inferred that 
the radio (whether from a jet or flaring events) could contribute 
significantly to the NIR data. However, emission in excess to that 
extrapolation was also seen in the optical bands, as suggested by 
the simple power law plotted through the radio data in Fig. [TT] 
This emission probably originates in both the disc and the com- 
panion star. Finally, a break was necessary to account for the 
X-ray/soft y-ray data, well fitted by a steeper slope than in ra- 
dio/NIR/optical, with a very soft power law photon index value 
(-1-2.34): the hot medium was thus detected. Therefore at least 
three components were needed to fit our broad band spectra, 
from radio to X/y-rays, probably accounting for the presence of 
matter from a previous ejection event, a disc (plus the companion 
star) and a hot medium. However, we did not have enough data 
in the required range (0.1-2 keV and in the optical) to precisely 
constrain the SEDs: the disc contribution was too small. 

The shape of the SE D is similar to the one of the transien t 
LMXB XTE J1118-H480 dChatv et alj|200l IZurita et al.ll2006h : 
although in the LHS during its 2005 outburst, the latter work 
gave IR and optical slopes with different a values, softening 
from 0.49 at the peak, not far from our value but distinct, as 
XTE J1818-245 was in another spectral state, to 0.25 (decay). 
In the case of flie 2000 outburst, the SED of XTE Jl 118+480 
from radio to X-rays has been explained as a combination of 
synchrotron radiation from a jet and a truncated optically thick 
disc, whereas models assuming advection dominated accretion 
flows alone underestimated the optical and IR fluxes (Zurita et 
al. 2006 and references therein). In 2005, discrepancies observed 
between the optical and IR SEDs suggested that the IR was dom- 
inated by a different component (a jet?) whereas in the optical, 
the authors were targeting the disc evolution (as we show in this 
work for XTE J1818-245). Linear fits to optical SEDs have also 
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been performed for other BH XNe in outburst (lHvnesll2005l) : re- 
sults are relatively uniform and all SEDs exhibit quasi power law 
spectra, with a ranging from 0.5-1.5 (compatible with our spec- 
tral in dex of 1), all steeper than the canonical Sy oc v'^^*. In par- 
ticular, HynesetiD (I2OO2I) observed XTE J 1 859+226 during the 
decline from its 1999-2000 outburst: in the soft state, the slopes 
derived were always between 0.1 and 0.5. This is slightly differ- 
ent from our value found in the SIMS: we took it closer to the 
peak. The authors found that the UV/optical/X-ray data - when 
detectable - could be fitted with a simple black-body model of 
an accretion disc heated by internal viscosity and X-ray irradi- 
ation, but the inner radius could not be well constrained. They 
concluded that the flat-spectrum synchrotron emission may be 
important in the IR and optical in this source. However, they did 
not exclude the alternative explanation that the IR excess could 
come from the cool outer disc. In our case, although we do not 
have enough broad-band and time coverage with our data to con- 
strain models as in the previous works, the contribution of the ra- 
dio (under a jet or flares, see Sect. 14. 3t in synchrotron emission 
up to the NIR/optical could be important, but another compo- 
nent (for example from the cooling disc and/or the companion, 
or even an irradiated disc, see Sect. 14. 11 and 14. 2t is necessary to 
account for the NIR/optical excess observed. 



4.3. Radio constraints and links witli X/y-rays 

The VLA radio observations of XTE J18 18-245 started around 

6 days after the ASM peak of the X-ray outburst. While the 
presence of a possible major flare associated with this outburst 
was not covered by these data, two radio flares were detected 

7 and 15 days after the soft X-ray peak. The first one occurred 
around MJD 53603 when the source was still in the SIMS; no 
spectral information is available for this flare. A couple of days 
later, when the source was in the harder part of the SIMS, it 
displayed the behaviour shown in Fig. |9] reminiscent of the 
radio flares displayed by GRS 1915-Hl05 lMirabel et al.|[l998l) . 
The lack of strictly simultaneous multi-frequency light curves 
prevents us from deriving the spectral behaviour of the source. 



The VLBA upper limits of 1 mJy beam"' at both 2.3 and 
8.4 GHz were obtained with data taken during the second hour 
of this VLA run and lasting two more hours. These upper limits 
could indicate either that the source has faded significantly, as 
suggested by the VLA 15 GHz upper limit, or that the source 
has slightly expanded, lowering the flux of compact emission 
with respect to the total emission seen at VLA scales. The lack 
of simultaneous VLA observations covering the whole VLBA 
run prevents us from reaching a firm conclusion on this respect. 
The second radio flare took place around MJD 5361 1, just before 
a transition to the HSS (see HID in Fig. |4]i. The radio spectral 
index (a = -0.8 + 0.1) obtained during this observation is typ- 
ical of optically thin synchrotron radiation, probably as a result 
of freely expandin g plasma blobs previously ejected (see, e.g., 
iFender et al.ll2004) . This suggests that multiple ejection events 
took place in the outburst of XTE J1818-245. Interestingly, the 
source was also detected on MJD 53617, five days after the start 
of the HSS. The spectral index is very negative, suggesting that 
this might be residual emission from the second flare. However, 
the large error bars due to the faintness of the source prevent us 
from being sure about this last result. 



5. Summary and conclusions 

During its 2005 outburst, XTE Jl 8 18-245 showed Hght curves 
and spectral evolution of a LMXB in outburst, probably with 
a BH as the compact object. Spectral parameters were typical 
of soft-intermediate and high-soft states, showing the usual de- 
crease in disc temperature, increase in inner disc radius, and the 
disc flux decreasing as the high-energy flux becomes compar- 
atively stronger. However, the track in the HID was not com- 
pletely standard, and it was accompanied by at least two ra- 
dio flares in the soft intermediate state, without strong hard X- 
ray changes. Even if there was an evolution in the power law 
index and flux, we did not observ e the source switch to the 
LHS as we did fo r XTE J 1720-3 18 (Cado lle Bel etaLll2004l) or 
XTE J1817-330 jCadoUe Bel et al. 2008) . 

Radio flares have now been seen in a large num ber of sources 
while transiting from HIMS to SIMS (see, e.g., iFender et al.l 
,2004) . Several models attempt to characterize the expelled ma- 
terial, the simplest of which is adiabatically expanding ma- 
terial eje cted by the sy stern , detectable from its synchrotron 
emission dvan der LaanI Il966h . In this model, a lag of several 
tens of minutes is expected between the high and low radio 
radio frequencies. Suc h a delay has now been seen in sev- 
eral sources (see, e.g. lMirabeletal.lll998t IFender et aP l2004t 
[Rodriguez et"ani2008allbl) . In the case of XTE J1818-245, the 
two radio flares were clearly associated with two different ejec- 
tion events. The fact we detected radio emission 5 days after 
the source entered the HSS is intriguing and can not be due to 
the residuals of the flare. One could compare this beha viour to 
XTE J 1650-500 or XTE J 1550-564 dCorbel et al.ll2004 . where 
the authors surprisingly detected a faint level of radio emission 
in the soft state. It was consistent with the emission of previ- 
ously ejected material interacting with the interstellar medium: 
this might also be the case here. 

In most BH binaries, a major ejection event occurs when 
the compact jet, present in the LHS in the very first stage of 
an outburst, turns off as the source enters the HSS. But in some 
cases the source then alternates between softer and harder states, 
and repeated ejections on a timescale of days to months may 
occur. This was the c ase for, e.g., A0620-00, GRO J1655-40 
and XTE J1859+226 (iHarmon et al.| [T99llKuulker s et al.lll999l: 
iBrocksopp et al.l2002l) . However, contrary to what was observed 
in these sources, we did not observe a clear correlation between 
the (hard) X-ray spectral behaviour and the occurrences of the 
radio flares in XTE J1818-245. Although some "switches" back 
and forth between SIMS and HSS were seen in the HID, this 
source remained in rather soft states, and only the last radio flare 
event was followed by a spectral state transition to a HSS, as in 
XTE J 1859+226 (Brockso pp et al. 2002). 

Our analysis clearly shows the importance of simultaneous 
multiwavelength and monitoring campaign of sources in out- 
burst. In particular, and although XTE J18 18-245 shows quite 
common behaviour for a BH in outburst, some specific points 
are quite intriguing. Besides, constraints from the modeling of 
the SEDs revealed that its behaviour was not fully compati- 
ble with the interpretation formed for other sources in terms of 
radio/X-ray evolution. The radio extrapolation might be well in- 
terpreted as non-thermal emission from ejected material, such 
as discrete ejection events that rapidly faded, except for resid- 
ual emission seen a few days after entering the HSS. The radio, 
whether emitted by a jet or as flaring events, could contribute sig- 
nificantly to the NIR data, but an excess was seen and could be 
explained by contributions from the cooling disc plus the com- 
panion. However, alternative processes, such as X-ray irradia- 
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lion, are sometimes inferred, depending on the source and data 
available, but there is no consensus yet on that issue (even for 
the same source), nor on the exact disc geometry /hard radiation 
origin. Clearly, the diversity and complexity of the interplay of 
the different emitting media during state transitions need more 
campaigns and observations to be better understood. 
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